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The exclusive electrodisintegration of the deuteron is considered within the Bethe-Salpeter 



^ ■ approach with a separable interaction kernel. The relativistic kernel of nucleon-nucleon 

^ interaction is obtained considering the phase shifts in the elastic neutron-proton scattering 

and properties of the deuteron. The differential cross section is calculated within the 
^ I impulse approximation under several kinematic conditions of the Bonn experiment. Final 

I state interactions between the outgoing nucleons are taken into account. Partial- wave 

^ ■ states of the neutron-proton pair with total angular momentum J = 0, 1 are considered. 

cn 

^ 1 Introduction 
(N 

The exclusive electrodisintegration of the deuteron is a useful instrument which makes it possible 
to investigate the electromagnetic structure not only of the bound state - deuteron but also the 
^ ' scattering states of the neutron-proton (np) system. Many approaches have been elaborated to 
■ describe this reaction for the last 40 years [1, 2, 3, 4, 5]. The simplest of them considered the 
electrodisintegration within a nonrelativistic model of the nucleon-nucleon (NN) interaction and 
outgoing nucleons were supposed to be free [1] (the plane- wave approximation - PWA). Those 
approaches were in a good agreement with experimental data at low energies. However, further 
investigations have shown that the final state interaction (FSI) between outgoing nucleons, two- 
body currents and other effects should be taken into account to obtain a reasonable agreement 
with existing experimental data at higher energies. Most of these effects have been considered 
within nonrelativistic models [2, 3]. In relativistic models, FSI effects could be calculated within 
quasipotential approaches using the on-mass-shell nucleon-nucleon T matrix [4, 5]. 

One of the fundamental approaches used to describe of the np system is based on the Bethe- 
Salpeter (BS) equation [6]. The separable ansatz [7] for the NN interaction kernel was success- 
fully applied to solve this equation. However, the approach was not applicable for the most 
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of high-energy NN processes for a while since calculated expressions contained nonintcgrable 
singularities. The problem was solved in [8, 9, 10] where kernels of special type were proposed. 
Using them FSI can be taken into account (in particular, considering the electrodisintegration) 
in a wide range of energy. 

In the present paper, the electrodisintegration cross section is calculated under different 
kinematic conditions of Bonn experiments [11, 12]. The rank-six NN interaction potential MY6 
[9] is used to describe scattered states ^Si-^Di and the deuteron. The uncoupled partial-wave 
states with total angular momentum J — 0,1 {^Sq, ^-Pi, ^Pq, ^Pi) are described by multirank 
separable potentials [10]. 

The paper is organized as follows. In Sec. 2, the exclusive three-differential cross section of 
the d{e, e'p)n reaction is calculated within the relativistic impulse approximation. The used BS 
formalism is presented in Sec. 3. The details of calculations are considered in Sec.4. Then the 
obtained relativistic results are compared with experimental data of the Bonn experiments [11, 
12] in Sec.5. 



2 Cross section 

When all particles are unpolarized the exclusive electrodisintegration of the deuteron d{e, e'p)n 
can be described by the differential cross section in the deuteron rest frame - laboratory system 
(LS), which has the following form: 

dE'jn'jnp ~ 8Mrf(27r)3 ./TTv\pp\ - Ep^cosOp ^ ' 

X Ko^oo + + W__) + 2ll_ cos 20 ReW+_ - 2/°_ sin 20 lmW+_ 

-2/g+ cos Rc{Wo+ -Wo-)- 2Z°+ sin lm{Wo+ + l^o-)] • 

where (TMott = (a cos |/2£^eSin^ |)^ is the Mott cross section, a = e^/(47r) is the fine structure 
constant; is the mass of the deuteron; q = Pe — p'e = ('^^ q) is the momentum transfer; 
Pe — {Eg, I) and p'^ — {E'^, l) are initial and final electron momenta, respectively; Jig is the 
outgoing electron solid angle; 9 is the electron scattering angle. The outgoing proton is described 
by momentum {Ep = yjp^^'^rrn^ , m is the mass of the nucleon) and solid angle Qp = {6p, (p) 
where 6p is the zenithal angle between the Pp and q momenta and is the azimuthal angle 
between the (ee') and (qp) planes. Factor rj = /s can be calculated through the np pair 
total momentum P squared: 

s = P^ = {Pp + Pnf = Ml + 2MdUJ + q\ (2) 

defined by the sum of the proton pp and neutron p„ momenta. The photon density matrix 
elements have the following form: 

where — —q^ is introduced for convenience. The hadron density matrix elements have the 
following form: 

W^x' = W^^e^^el,, (4) 
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where A, A' are photon hehcity components [13] , can be calculated using the photon polarization 
vectors e and Cartesian components of hadron tensor 

W^^-^IYI \< '■ ^^s\jM : IM >|', (5) 

where S is the spin of the np pair and Ms is its projection; s^, s„ and Sp are deuteron, neutron 
and proton momentum projections, respectively. The matrix element < np : SMs\3ij\d : IM > 
can be written using the Mandelstam technique [14] as follows: 

CM 

-^Sp {A(£-^)V;5M,(P°",P*'^"; X (6) 

r(n)(^)^(n) _ (-1)>- I) r- (p+ (-l)"|;i^(o))} 

within the relativistic impulse approximation in LS. The sum over n = 1,2 corresponds to the 
interaction of the virtual photon with proton and with neutron in the deuteron, respectively. 
The total P'^^ and relative p*^'^ momenta of the outgoing nucleons and the integration mo- 
mentum p^^ are considered in the final np pair rest frame - centcr-of-mass system (CMS), p 
denotes the relative np pair momentum in LS. To perform the integration, momenta p, q and 
deuteron total momentum i^(o) — {Ma, 0) in LS are written in CMS using the Lorenz-boost 
transformation C along the q direction. The np pair wave function ipsMs is transformed from 
CMS to LS applying the corresponding boost operator A. A detailed description of i^sMs, the 
nth nucleon interaction vertex F^"-*, the propagator of the nth nucleon S^'^\ and the deuteron 
vertex function can be found in our previous works [15, 16]. 

3 Separable kernel of NN interaction 

The outgoing np pair is described by the T matrix which can be found solving the inhomoge- 
neous Bethe-Salpeter equation [6]: 

T{p\p-P)^V{p\p-P) + ^ j d'kV{p',k;P)S2{k;P)T{k,p;P), (7) 

where V is the NN interaction kernel, 5*2 is the free two-particle Green function: 

S^\k; P)={lf + ]^-mf\\f-]/,-m) (8) 

and p {p') is the relative momentum of initial (final) nucleons, P is the total np pair momentum. 
To solve the BS equation (7) partial-wave decomposition [17] for the T matrix: 

Tap,'rsip',P; Pio)) = Yl ^'^M, |p'|;Po, \p\\s){yaM{-p')Uc)ap ® iUcylMiP))s'r (9) 

JMab 

is used. Here P(o) = {y/s, 0) is the np pair total momentum in CMS, Uc = ^7^7° is the charge 
conjugation matrix. Indices a, b correspond to the set ^^'^^Lj of spin S, orbital L and total J 
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angular momenta, p = + defines a positive-energy partial-wave state, p = — corresponds to a 
negative-energy one. Greek letters {a,/?, 7, 5} in (9) are used to denote Dirac matrix indices. 
The spin-angle functions: 



mLmsmim2pip2 



arc constructed using free nucleon Dirac spinors u, v. It should be mentioned that only positive- 
energy states with p — + are considered in this paper. Performing similar decomposition for 
V, the BS equation for radial parts of the T matrix and kernel V is obtained: 

tab{p'o,\p'\;PoAp\;s) =Vab{p'o,\p'\;po,\p\;s) (10) 
+00 00 



^ I J^'^(^\^\'"o.cip'o,\p'\;ko,\k\;s)Scdiko,\k\;s) tdb{ko,\k\;po,\p\;s). 



-00 



47^3 

To solve the resulting equation (10), a separable ansatz [7] for the interaction kernel v is used: 

N 

Vabip'o, |p'|;Po, \p\;s) = Xij{s)g^^\p'o, \p'\)gf\po, \p\), (H) 



where is a rank of a separable kernel, Qi are model functions, is a parameter matrix. 
Substituting v (11) into BS equation (7), we obtain the t matrix in a similar separable form: 



N 

Ml 



where: 



and 



tabip'o: |P'|;P0, \P\;S) = Tij{s)gP{po, \p'\)gf\pO: \p\) (12) 

«J=1 



njis)^l/{\-^is) + h,,{s)), (13) 



(14) 



are auxiliary functions, = \/k'^ + w?. Thus, the problem of solving the initial integral BS 
equation (7) turns out to finding functions gi and parameters Ay of separable representation 
(11). They can be obtained from a description of observables in np elastic scattering [9, 10, 18, 
19, 20]. 

4 Final state interaction 

To calculate FSI we need to consider the interacting np pair. The outgoing nucleons are 
described by the BS amplitude which can be written as a sum of two terms: 

V'5m,(p,P*;P) = 4m,(p,P*;^) + 4^'52(p;P) j d'kV{p,k;P)i^sMs{k,p*;P). (15) 
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The first term 

4m«(p,p*;^) = (27r)Ws(p;^)5(p-p*) (16) 

is related to the outgoing pair of free nucleons, XsMs ^ spinor function for two fermions. The 
second term in (15) corresponds to the final state interaction of the outgoing nucleons. If we 
use the following relation: 

d'k V{p, k- P)^sMs{k,P*; P)^ j d'k T{p, k; P)4ts(^>P*; P) 
then it can be transformed into 

41^(p,p*;P) = 4mS2{p;P)T{p,p*;P)xsMsiP*;P), (17) 

here (t) means that this part of the np pair wave function is related to the T matrix. Using 
the partial- wave decomposition (9) expression (17) can be written as follows: 

^^sls(P^P*-,P)-^^^ E Ci^SMsyUPlyaM{p)cl>a,J:Ls4P0:\p\;s), (18) 

LmJMa 



where p* = {0,p*) with = y s/4 — is the relative momentum of on-mass-shell nucleons 
in CMS, p* denotes the azimuthal angle 6p* between the p* and q vectors and zenithal angle 
(p. Since only positive-energy partial-wave states are considered here the radial part is: 

, ^ I I ^ ta,J:LS+{PO,\p\;0,\P*\;s) 

<Pa,J:LS^ipo, H s) = (^/2 _ + _ ^2 " (19) 

According to definition (10) spin-angle functions y can be written as a product of Dirac matrices 
7 in the matrix representation [15] as follows: 

matrices QaM are shown in Table 1. Decomposition (18) is considered in detail in [22]. Using 
definition (15) and substituting (16), (18) into (6), the final expression for hadron current 
< np : SMs\jix\d : IM > can be obtained. It consists of two parts. One of them: 

< np : SMslJM : IM >W= t J] {A(£-^)x5m, (p*""; P^ A(>C) x (21) 

n=l,2 

r(n)(^)^(n) (^K^ _ _ |^ (^^. + (_l)n|. ^^^^ j | 

corresponds to the electrodisintegration in PWA. Another one: 

< np : SMs\j„\d : IM (22) 

oo oo 1 27r 

4^ E E ci^sMsyuPi I dpr I {p'^n'd\p--\ I dcos^- / dct> x 

n=l,2LmJMjL'lm' _^ q _i q 

Sp{A(£-^)j),^,5M,(p^-)A(£)r(r)(g)x S^-^ - (-1)>- |) yiiSm' (p+ (-1)"|)} 

n'L(Po",|p°"l;0, \p*\;s) 



X 



{^/2-E, + ie)^-pl 
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^^/(po + (-i)"|,P+(-i)"f;i^(o)) 







QaM 


3 p+ 




-75 

■^m)75 

-2^(1^1-^2) 


3p+ 


fs 1 

V 2|PI 


ij>l<M) - Uuifi 


3L>+ 


73 U + lp^iPi-^M){p, 



Table 1: Spin-angular parts QaM (20) for the np pair; pi = {Ep,p), p2 = {Ep, —p) are on-mass-shell 
momenta, Ep = y^p^"!]!^^. ^ matrices are defined as in [21]. 



corresponds to the process when FSI is taken into account. Here gi denotes the radial part 
of the deuteron vertex function T'^ . The part Sp{. . .} has been calculated using the algebra 
manipulation package MAPLE. The three-dimensional integration over p^^, \p^^^\ and cos^^"^ 
has been performed numerically using the programming language FORTRAN. 



5 Results and discussion 

The calculations of the deuteron electrodisintegration within PWA were considered in [9, 16]. 
As it was shown in [3, 23], a contribution of FSI effects increases with increasing nucleon 
energies or/and momentum transfer. The reaction near the threshold was considered in [23] 
(Sacley experiments, see [24, 25]). 

In the paper the differential cross section (1) has been calculated under five kinematic 
conditions of the Bonn experiments [11, 12] (described in Table 2) and is shown in Figs. 1-5. 
The calculations have been performed within the relativistic impulse approximation for two 
different cases: when the outgoing nucleons are supposed to be free (PWA) and when the final 
state interaction between the nucleons is taken into account (FSI). The partial- wave states of 
the np pair with total angular momentum J — 0,1 have been considered. The used relativistic 
model consists of two parts: the separable potential MY6 [9] for the bound (deuteron) and 
scattered '^Si-^Di states and separable potentials of various ranks [10] - for all other partial- 
wave states (^5*0, ^Pi, ^Pq, ^Pi). 

In Figs. 1-5, relativistic PWA (sohd red line) and FSI (dashed blue line) calculations are 
shown as a function of ^/s. 

As it is seen from Figs. 1-5, the effect of FSI increases the cross section for kinematic condi- 
tions [11, 12]. 

In Fig. 2, the calculation with FSI describes the experimental data on ^/s from 2.15 GeV to 
2.3 GeV and differs from experimental data starting from 2.3 GeV. Calculations of the cross 
section shown in Figs. 1-2 require additional investigations at high s (for instance, the influence 
of non- nucleon degrees of freedom). 

In Figs. 3-5, PWA and FSI calculations under the kinematic conditions [12] are shown. The 
effect of FSI is small at the threshold but it becomes bigger with increasing ^/s. 
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The inelasticity effects are nonzero for kinematic conditions of Bonn experiments [11, 12]. 
However, thev will be discussed in a senarate Daner. 




Figure 1: Cross section (1) depending on ^Figure 2: The same as in Fig.l but under kine- 
calculated under kinematic conditions set I of thematic conditions set II of the Bonn experiment 
Bonn experiment [11]. i^'^]- 




Figure 3: Cross section (1) depending y/s cal- Figure 4: The same as in Fig.l but under kine- 
culated under kinematic conditions set III of thematic conditions set IV of the Bonn experiment 
Bonn experiment [12]. [12]- 
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Vs, GeV 

Figure 5: Cross section (1) depending on ^/s 
calculated under kinematic conditions set V of 
the Bonn experiment [12]. 
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Table 2: Kinematic conditions considered in the paper. Here all quantities are given in LS. In addition 
to those which are defined in the text, they are: angle 6qe between the beam and the virtual photon; 
neutron momentum p„ and angle 9n between the neutron and the virtual photon (pp, Op - the same 
for the proton); 6pe {6qe) - the angle between the beam and the proton (virtual photon). 
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